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ABSTRACT 

We use a cosmological galactic evolutionary approach to model the Milky Way. 
A detailed treatment of the mass aggregation and dynamical history of the growing 
dark halo is included, together with a self consistent physical treatment for the star 
formation processes within the growing galactic disc. This allows us to calculate the 
temporal evolution of star and gas surface densities at all galactic radii, in particular, 
the star formation history (SFH) at the solar radius. A large range of cosmological mass 
aggregation histories (MAHs) is capable of producing a galaxy with the present day 
properties of the Milky Way. The resulting SFHs for the solar neighbourhood bracket 
the available observational data for this feature, the most probable MAH yielding the 
optimal comparison with these observations. We also find that the rotation curve for 
our Galaxy implies the presence of a constant density core in its dark matter halo. 

Key words: Cosmology: dark matter — Galaxy: evolution — Galaxy: halo — galax- 
ies: formation — solar neighborhood — stars: formation 



1 INTRODUCTION 

, The star formation history (SFH) of galaxies is a key ingre- 
' dient in describing their evolution. From an observational 
point of view, the SFH can be studied by measuring the 
1 integrated light of the galaxy population as a whole or of 
' a homogeneous sample of galaxies at different redshifts, at 
wavelengths associated with massive star formation (e.g., 
Madau, Pozzetti & Dickinson 1998; Lilly et al. 1998). An 
alternative method, which avoids the difficulties in selecting 
a homogeneous sample of galaxies, determines the SFH of a 
specific galaxy using multi-band photometry and chemical 
abundances combined with population synthesis models. In 
the case of the Milky Way (MW), the large amount of avail- 
able data aUows to qualitatively constrict the input SFHs of 
these methods with some accuracy (e.g., Prantzos & Aubert 
1995; Chiappini, Matteucci, & Gratton 1997; Carigi 1997; 
Boissier & Prantzos 1999). 

On the same line of multi-band photometry, but with 
higher precision, the SFH of a system can be obtained 
through the study of its resolved stellar population. This al- 
lows to construct colour-magnitude diagrams, which contain 
detailed information on the SFH behind such diagrams. In- 



verting the problem however, is not trivial, with the answer 
depending sensitively on the assumed metallicity of the ob- 
ject being studied at every time, i.e., the enrichment history. 
The accuracy of these inferences being also crucially deter- 
mined by the depth and error level of the observations, these 
requirements have so far severely constrained the applicabil- 
ity of such direct inferences to only a handful of astrophysical 
systems. Still, such methods have recently been applied to 
systems where the metallicity has been independently deter- 
mined, and where high quality data exist. Examples being 
the works of Chiosi et al. (1989), Aparicio et al. (1990) and 
Mould et al. (1997) using Magellanic and local star clusters, 
and Mighell & Butcher (1992), Smecker-Hane et al. (1994), 
Tolstoy & Saha (1996), Aparicio & Gallart (1995), Mighell 
(1997) and Hernandez et al. (2000a) using local dSph's. 

With the coming of high quality photometric data from 
the Hipparcos satellite, colour magnitude diagram inversion 
methods using rigorous statistical analysis (Tolstoy & Saha 
1996; Hernandez et al. 1999 and more references therein) can 
be applied to the solar neighbourhood SFH. For example, 
the advanced Bayesian analysis introduced by Hernandez et 
al. (2000b), allows the recovery of the underlying SFH with- 
out the need of assuming any a priori structure or condition 
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on the SFH. With the Hipparcos catalogue the time resolu- 
tion of this technique is ^ 0.05 Gyr, however, the age range 
which it allows to explore is small (the last 3 Gyr); this 
limitation is related to the completeness of the Hipparcos 
solar neighbourhood sample. With a standard parametric 
maximization technique, Bertelli & Nasi (2001) were able to 
recover the SFH along the whole life of the solar neighbour- 
hood, at the expense of a low time resolution. 

The high quality SFH which can be inferred locally al- 
lows the unique opportunity of studying directly the evolu- 
tion of an individual cosmological object (the Galaxy). This 
avoids the problems inherent to high redshift studies, such 
as the uncertainties in relating a population of poorly un- 
derstood and barely resolved objects to a particular class of 
present day systems. 

Once the SFH of an observed system is known, the goal 
is to understand its physical evolution. From a theoretical 
point of view, the challenge is to link the SFH to the struc- 
ture, dynamics and hydrodynamics of modeled galaxies. Ex- 
amples of galaxy evolutionary studies in a cosmological con- 
text are: White & Frenk (1991), Lacey & Silk (1991), Kauff- 
mann. White & Guiderdoni 1993, Cole et al. (1994), Baugh, 
Cole & Frenk (1997), van don Bosch (1998), Somorvillc & 
Primack (1999), Buchalter, Jimenez & Kamionkowski (2001) 
and Kauffmann, Chariot & Balogh (2001). A potential short- 
coming of many of the above approaches remains the use of 
somewhat ad hoc or empirical schemes for calculating the 
star formation (SF) in the forming discs. 

We improve on this last point by including a physi- 
cally self consistent scheme of calculating the SF process, 
where the SF in the disc is induced by gravitational insta- 
bilities and is self-regulated by an energy balance in the 
ISM (Firmani & Tutukov 1992, 1994; Firmani, Hernandez, & 
Gallagher 1996). Combining the above with a cosmological 
evolutionary approach (Avila-Reese, Firmani & Hernandez 
1998 and Firmani & Avila- Reese 2000) yields a powerful tool 
for exploring the consequences of any particular cosmology 
in terms of a detailed SFH (Avilar Reese & Firmani 2001). 

In this work we make use of this evolutionary approach 
to infer the SFH of the MW. The cosmological scenario used 
allows to constrain the SF history of the solar neighbour- 
hood to a narrow range of possibilities, and to identify the 
most likely of them. Comparison with the available direct 
observational inferences shows that the local SFH has in 
fact been very close to the maximum likelihood solution we 
find. Along similar lines, Hernandez & Ferrara (2001) have 
studied the link between the cosmological build up of the 
Milky Way and the IMF of population HI stars, by anar 
lyzing the mctallicity distribution of extremely metal poor 
Galactic halo stars. 

In §2 a brief description of the method is presented, 
emphasizing on the physics of the SF. In §3 we apply this 
method to model a MW galaxy. The SFH at the solar radius 
is presented and compared with observations in §4. Section 
4.1 is devoted to exploring the sensitivity of the predicted 
SFH to the input parameters which define our MW model; 
these parameters were settled from observational constraints 
which are subject to uncertainties. In the discussion we com- 
pare our results with those of chemo-spectrophotometric 



models (§5.1), and wo discuss the possibility of an inter- 
mittent SF, as well as other physical mechanisms capable 
of triggering SF (§5.2). Our conclusions are presented in §6. 
Throughout this work we assume the currently favored cos- 
mology, a flat universe with non-vanishing vacuum energy 
{Qm = 0.3, = 0.7, Ho = 65kms-^Mpc-^,a8 = 0.9). 



2 THE MODEL 

The overall method of disc galaxy evolution used here was 
presented in Firmani & Avila-Reese (2000, henceforth FAOO) 
and Avila-Rccsc & Firmani (2000, henceforth AFOO) (see 
also Firmani, Avila-Reese & Hernandez 1997). Galactic discs 
form inside-out within a growing CDM halo, as the baryonic 
content of the growing halo is incorporated into the disc. 
The evolution of the halo is set by its mass aggregation 
history (MAH), determined by an extended Press-Schechter 
formalism. For a given present-day total virial mass Mv, a 
statistical set of MAHs is calculated from the initial Gaus- 
sian density fluctuation field once the cosmology and power 
spectrum are fixed (Lacey & Cole 1993; Avila-Reese et al. 
1998). The behaviour of the growth of the halo in mass as 
a function of redshift can be seen in Fig. lb, which shows a 
sample of five realizations of possible MAHs. Fig. la gives 
the average over a statistical sample of 2 x 10* Monte Carlo 
realizations for a 2.8 x 10^^ Mq halo. Although the "average 
MAH" clearly smoothes over the discontinuous behaviour of 
any particular MAH, it provides the simplified version for 
the most probable mass aggregation process of a large sam- 
ple of galaxies. Further, the discontinuous behaviour seen in 
the individual MAHs of Fig. lb, as refiected in gas accretion, 
will be naturally smoothed over by the gas virialization and 
cooling processes in the growing galactic halo. 

The density profile of the virialized part of the grow- 
ing halo is calculated with a generalization of the secondary 
infall model, based on spherical symmetry and adiabatic in- 
variance, but allowing for non-radial orbits with fixed peri- 
centre to apocentre distance ratios e (Avila-Reese et al. 
1998). The halo density profile we obtain depends on the 
MAH and e, with the average MAH and e « 0.2 giving den- 
sity profiles very close to the Navarro, Frenk & White (1997, 
NFW) profile (Avila-Reese et al. 1998, 1999). 

A fraction fa of the mass of each collapsing spheri- 
cal shell is transferred in a virialization time into a central 
disc gas layer. Assuming angular momentum conservation, a 
given gas element of the shell falls into the equatorial plane 
at the position where it reaches centrifugal equilibrium. Each 
shell is assumed to be initially in solid body rotation, with 
a specific angular momentum given by: 

GM.(^)^/V 5 1 d\E{t)\ \ 

^'"^ ^ \E{t)\'/' V2Mv(t) 2\E{t)\dM^{t)J' 

where E is the energy of the halo at time t, and A is the 
halo spin parameter assumed here to be constant in time. 
The contraction of the dark halo due to the gravitational 
drag of the infalling gas is calculated within am adiabatic 
invariance scheme. A nearly exponential disc surface density 
arises naturally in this inside-out scheme of disc formation 
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Figure 1. a) Average MAH for a galaxy having the present day 
asymptotic circular velocity, V50, as the Milky Way. b) Five ran- 
dom realizations of cosmological MAHs yielding the same final 
asymptotic circulaj: velocity as in a). 

(e.g. FAOO). The A parameter determines the scale length 
and surface density of the disc. According to analytical and 
numerical studies, A has a log-normal distribution with an 
average of ~ 0.05 and a width in the logarithm of less than 
1.0 (e.g., Catelan & Theuns 1996 and references therein). 

Wc note that if the dynamical consequences of the 
MAIfs we are using were taken into account, the repeated 
merging of substructures would destroy the forming discs 
(Navarro & Steinmitz 2000). This is a shortcoming gen- 
eral to all hierarchical galactic formation models, related 
to the overabundance of substructure seen at galactic scales 
(Klypin et al. 1999, Moore et al. 1999a), the excessive mass 
concentration in the central region of CDM haloes (e.g., 
Moore 1994, Burkert & Silk 1997, Hernandez & Gilmore 
1998), and the loss of angular momentum of the gas dur- 
ing incorporation into the forming galaxy (Navarro & Bonz 
1991). All these difficulties might indeed imply modifications 
to the overall standard scenario. At this point we introduce 
as a working hypothesis the assumption that such modificar 
tions might result only in a more gentle and homogeneous 
mass accretion, without altering the overall structure for- 



mation and growth model, which appears at the moment 
a very robust prediction of the Gaussian initial fluctuation 
hypothesis within the inflationary CDM cosmology. 



2.1 Disc star formation 

As the disc begins to form in the centre of the growing CDM 
halo, transformation of cold gas into stars takes place when- 
ever Toomre's instability parameter for the gaseous disc, 
QsC^) = f's(^)'t(^)/^g(^)5 is less than a given threshold, q 
(where Vg and Eg are the gas turbulent rms velocity and 
surface density, and k is the epicyclic frequency, Toomre 
1964). A fraction of the formed stars explodes as SNae, re- 
heating the gas disc column and increasing Qg{r) above the 
instability threshold. This inhibits SF temporarily, whilst 
the gas disc column dissipates its turbulent energy, decreas- 
ing Vg back to the value when Qgir) becomes smaller than 
q. Actually we calculate the SF rate (SFR) directly from the 
equation that relates the energy input with the energy dissi- 
pation assuming that Qg (r) is always equal to q at all radii, 
i.e. we allow only for the stationary solution (a more de- 
tailed description can be found in Firmani & Tutukov 1994; 
Firmani et al. 1996). We also consider heating due to the 
infalling gas, although its contribution is much smaller than 
the SN input. 

While the SFR is rather insensitive to q, the gas disc 
height strongly depends upon it. Numerical and observa- 
tional studies suggest thresholds of the order of g = 2 (see 
references in FAOO). Further, for g = 2, we obtain realistic 
values for the gas disc height of a Milky Way (MW) model 
at the solar neighbourhood. The turbulent dissipation time 
td{r) is calculated following a simple approach given in Fir- 
mani et al. (1996): 



td 



6r 

a — = a 



1 



2n 



V/r + dV/dr 



(2) 



where 5r is the mean free radial distance traveled by a tur- 
bulent cell at Galactocentric radius r with a speed Vg {5r is 
assumed to be much smaller than the orbital length) , V and 
Q. are the circular and angular velocities, and a is a param- 
eter close to unity. Estimates of the turbulent dissipation 
time in the solar neighbourhood obtained in compressible 
magneto-hydrodynamic simulations of the ISM (Avila-Reese 
& Vazquez-Semadeni 2001) give a lower limit that is roughly 
half the time given by Firmani et al. (i.e. a > 0.5 in cq. (2)). 
It should also be taken into account that in the real ISM, be- 
sides turbulence, magnetic fields and cosmic rays contribute 
to the pressure, in such a way that the td found in the tur- 
bulence simulations could be only a lower limit to the global 
energetic dissipation time in the ISM. We shall use here 
Q = 0.7 

The evolution of the stellar populations is followed with 
a parameterization of simple population synthesis models 
(sec Firmani & Tutukov 1992, 1994). A Salpotor IMF with 
minimal and ma^ximal masses of O.IM© and lOOM©, and 
solar metallicities were used. 

The aaimuthally averaged dynamics of the evolving gas 
and stellar discs coupled with the dark halo are treated by 
solving the corresponding hydrodynamical equations. With 
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our method we obtain the temporal evolution of the gas 
infall rate on the disc Eg(r, i), the SFR Es(j", t), the disc 
gas and star surface densities, Eg(r, and Es(r, t), and the 
circular velocities due to the disc and halo, Vd(r, and 
Vh(r, t), respectively at all radii. In our model the halo struc- 
ture and the disc gas infall rate are directly related to the 
cosmology. 

We find that the local SFR typically follows a Schmidt 

law, Es('') oc Sg", with n « 2 along a major portion of the 
disc (Firmani ot al. 1996; AFOO). This SF law is basically 
a consequence of the self-regulation mechanism we use. It 
is important to mention that in our models, self-regulation 
takes place only within the disc and not at the level of a 
hypothetical intra-halo medium in virial equilibrium with 
the CDM halo (e.g.. White & Frenk 1991; Cole et al. 1994; 
Baugh et al. 1997; van den Bosch 1998; Somerville & Pri- 
mack 1999). The ISM of normal disc galaxies is dense and 
very dissipative. Hence, one expects gas and kinetic energy 
outflows to be confined mostly near the disc (see numeri- 
cal results and more references in Avila-Reese & Vazquez- 
Semadeni 2001). Nonetheless, in some cases and at early 
epochs, the disc-halo feedback may be important; in our 
model this ingredient may be taken into account by the free 
parameter f^, at least in a first approximation. 



3 MODELING THE MILKY WAY 

For a given cosmology a model is completely determined by 
the virial mass Mv, the MAH, the spin parameter A, and 
the disc mass fraction fj, as all other parameters are fixed 
by physical considerations. 

In order to calculate a model representative of the MW, 
we have to correctly chose these input parameters so as to 
fit the observational constraints. Fortunately, each one of 
these parameters is tightly related to a given present-day 
MW feature, although with some (small) degeneracy. Due 
to this degeneracy and to the observational uncertainties in 
the MW observational constraints, it is necessary to carry 
out some "parameter-space" exploration. 

The total mass of the MW is amongst the most poorly 
known of all Galactic parameters. Several estimates give val- 
ues around 1.0 — 4.0 x lO^^M© (see e.g., Wilkinson & Evans 
1999 and more references therein). A more robust quantity 
is the mass within a certain large radius (50 kpc for exam- 
ple) constrained by the motions of satellite galaxies, globular 
clusters, the local escape velocity of stars, etc. Using these 
constraints together with the observed rotation curve of the 
disc, Kochanek (1996) determined the MW mass within 50 
kpc. This mass implies a circular velocity V50 of 2061J5 
km s"i(206l25 km s"\ at 90% C.L.). A more recent deter- 
mination of V50 by Wilkinson & Evans (1999), constructed 
from the current data set of objects (27) with known dis- 
tances and radial velocities at Galactocentric radii greater 
than 20 kpc, is consistent with the values found by Kochanek 
(1996). A V50 of 206lJ° km s~^is the observational con- 
strain we use to fix the virial mass of the MW model for 
a given MAH. Since the MAH determines the halo density 



Table 1. Properties of the MW: observations and model results 



Observable Predicted value " Observed value Reference 



Vc(50) [km s-i] 


208 


206 ± 10 


1,2 


rs[kpc] 


3.0 


3.0 ± 0.5 


3 


Ms = [IQio Mq] 


4.4 


4.5 


4,5 


Mv [10^2 Mq] 


2.8 


1-4 


2,6 


Vmax'* [km S-1] 


235 


220 ± 10 


7,8,9,5 


Lb [IQIOLsq] 


1.7 


1.8 ± 0.3 


10,9 


r_g [kpc] 


4.3 


4-5 


11 


So,K ^ [Lx0PC-2] 


810 


1000 ± 200 


11 


fg 


0.23 


0.15-0.20 


12 


SFR [Moyr-i] 


2.9 


2-6 


13,14 


Solar neighborhood 


Es[M0pc-2] 


41.3 


43 ± 5 9 


4 


Sg[M0pc-2] 


13.6 


13 ± 3 


15 


E4Mo,Gyr-i pc-2] 


3.1 


3-5 


16,17 


LI - K [iiiixg] 


:-!.l5 




8 



" The MW model was obtained tuning three input parameters 
in order to reproduce the first three quantities, which are 
constraints and not predictions 
' Circular (asymptotic) velocity at 50 kpc. 
Stellar (disc+bulgc) mass. 

This maxiirmiri rotation velocity Vmax docs not take into 
account the nuclear and bulge region and it is for Galactocentric 
radii smaller than 15 kpc; the observational determinations at 
larger radii are uncertain. 
" B— band disc scale length. 

^ Disc central surface brightness in the JC— band, taJfing a mass 
to light ratio of 1.0 for this band. 

B This estimate includes the contribution of stellar remnants. 

References: 1. Kochanek 1996; 2. Wilkinson & Evans 1999; 3. 
Sackett 1997; 4. Mera et al. 1998; 5. Dehnen & Binney 1998a; 6. 
Mendez et al. 1999; 7. Fich & Tremaine 1991; 8. Binney & 
Merrifield 1998; 9. Binney & Tremaine 1987; 10. van der Kruit 

1986; 11. Kent et al. 1991; 12. sec BP99; 13. Prantzos & Aubert 
1995; 14. Pagel 1997; 15. Kulkarni & Heiles 1987; 16. Talbot 
1980 17. Rana 1987 



profile, Mv will be different for different MAHs; for the av- 
erage MAH, this mass is Mv = 2.8 X 10^2 Mq. 

The spin parameter A mainly influences the stellar disc 
scale length and the surface density. The stellar scale length 
rs of the MW inferred from observations in infrared bands 
or estimated from dynamical constraints is typically 3.0 ± 
0.5 kpc (Sackett 1997), this becomes our second constraint. 
However, see for example Drimmel & Spergel 2001 where 
values even smaller than those wore found. For the average 
MAH, the observed central value of ts yields A = 0.02. 

Once Mv is given, U= Ma/Mv is fixed by the MW 
disc-|-bulge mass Mj. Most baryons are in the disc-|-bulge 
system: 4 — 5 x 10^° M© in form of stars (Ms, our third 
constraint) and at least ~ 0.6 — 1.0 x 10^° Mq in form of 
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gas (Mg) (see Table 1). The estimated masses 
halo (« 4.6 X 10'' Mq, Haud & Einasto 1989) ar 
satellite galaxies (~ 6 — 9 x 10^ M©) represent 
contribution with respect to the disc+bulge sys 
the average MAH we obtain fd=0.021. 

The influence of different MAHs is felt upoi 
galaxy color index, gas fraction fg, and scatter 
Fisher relation. Unfortunately the MAH cannot 
by one-parameter and the color index and gas 
not accurately determined for the MW. Thus 
with a large range of possible MAHs, all of w 
models with the same V50, rs and Ms for the M 
the SFH in the disc of the evolving galaxy, and 
at the solar radius, is always highly sensitive 
MAH. The range of MAHs allowed by the sti 
straints taken hence defines a range of possil 
the solar neighbourhood, with the most likeb 
associated to the average MAH. A comparison 
ferred SFH of the solar neighbourhood, explo: 
fact shows the MAH of our Galaxy to have beei 
average one for systems of the MW mass. 

The observational estimates of V50, rs ai 
in Table 1 are thus the constraints we require 
a model to bo considered representative of the 
given MAH the previous conditions allow to de 
A and fd. When the average MAH is assume 
Mv= 2.8 X 10^^ A/0, A = 0.02 and fd=0.021. ' 
tions determine our MW central model, which 
cosmological maximum likelihood solution, givi 
tural parameters of our Galaxy (fiducial model 

It should be noted that the values of fd 
in order to reproduce the stellar mass and sc; 
the MW, are smaller than those used in other ! 
canton, Spcrgel & Summers 1997; Mo, Mao & 
FAOO, AFOO). The disc mass fraction in hak 
in these works was w 0.05, only slightly smal 
baryon to total mass fraction, Qb/^im, for the 
value of Qth^ ^ 0.010 (Hogan 1998) and Qm 
disc mass fraction increases to fb ~ 0.15 tal 
count the higher estimate of 0,bh^ («0.019) i: 
the deuterium absorption in quasars (Buries 
Thus, for the MW the baryon fraction (Md + il 
Msatciiitcs + ...)/Mv) seems to be smaller than uiiu umvei- 
sal baryon fraction. Cosmological numerical simulations in- 
deed show that most of the baryons in the past and in the 
present are not within gravitationally bound systems, but 
in form of diffuse and warm-hot intergalactic gas (Gen & 
Ostriker 1999; Davet al. 2000). Smaller values of fd allow for 
smaller values of A (given a fixed rs) that otherwise could 
lead the disc to destructive gravitational instabilities (see 
Mo et al. 1998; FAOO). If the distribution of A is mainly re- 
sponsible for the galaxy distribution in surface brightness, 
then, that high-surface brightness MW-typc galaxies should 
have A « 0.02 agrees with a distribution where low surface 
brightness (LSB) galaxies (A ^ 0.05) are equally abundant 
(e.g. McGaugh 1996). To conclude, the values for fd and A 
inferred for the average MAH seem to be consistent with 
several observational and theoretical pieces of evidence. 




R [kpc] 

Figure 2. a) Rotation curve decomposition for the fiducial 
Galaxy model giving the total, halo and disc contributions, solid 
dashed and dotted thick lines, respectively. The thin lines are 
analogous to the thick ones, showing the case of artificially im- 
posing a 10 kpc constant density core on the dark halo, b) Dif- 
ferent final rotation curves, for MW models resulting from the 
various MAHs shown in Fig. lb, with line style corresponding to 
the particular cases of Fig. lb. 



3.1 Rotation curve decomposition 

In Fig. 2a the rotation curve decomposition of the fiducial 
MW model is shown by the thick lines, with the continuous 
curve giving the total rotation curve, the dashed one the 
halo component, and the dotted one the disc contribution, 
at the present epoch. The maximum of the final rotation 
curve occurs at at « 2.2 is with a value of 235 km s~^. 

The model rotation curve is somewhat more peaJsed 
than the observed one (Fich & Tremaine 1991) and the halo 
component dominates over the disc component at almost 
all radii. This is a direct consequence of the high central 
mass concentration of the CDM halo. There is no way of 
obtaining a flat rotation curve for the parameters of the 
MW when a GDM halo is used. The density profile of our 
"virgin" CDM halo (without the gravitational drag due to 
disc formation) is well described by the NFW profile. Recent 
very high-resolution simulations show that the inner density 
profile of GDM haloes is even steeper than the slope of 
NFW profiles (Moore et al. 1999b; see also Fukushige & 
Makino 1997). The rotation curve of a MW model in this 
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T [Gyr] 

Figure 3. a) The thick solid curve gives the predicted SFH 
for the solar neighbourhood corresponding to the most probable 
cosmological MAH (the average MAH), with the corresponding 
gas accretion history at the solar neighbourhood given by the 
dot-dashed curve. Also shown are the observational inferences of 
BDBOO for the age of the solar neighbourhood, horizontal shaded 
box, and of Talbot 1980 and Rana 1987 for the present day value 
of the SFR, vertical shaded box. The jointed circles give the in- 
dicative inferences of Bertelli & Nasi (2001). Our ijrcdiction for 
the average MAH clearly complies with the available observa- 
tional restrictions for the local SFH. b) This panel shows the 
evolution of the gas and stars surface densities at the solar neigh- 
bourhood, dashed and solid lines, respectively. The present day 
values being in good accordance with observations. 

case would be even more peaked than that shown in Fig. 2, 
which should be taken as a lower limit in terms of the cen- 
tral overshooting. The intrinsic scatter in the density profiles 
of CDM haloes does not help to alleviate this problem, as 

can be seen in Fig. 2b, whore the total final rotation curves 
corresponding to the five particular MAHs of Fig. lb are 
shown. 

Several pieces of evidence point out that dark mat- 
ter haloes at all scales indeed have constant density cores 
(Moore 1994, Florcs & Primack 1994 and Burkcrt 1995 in 
dwarf galaxies, de Blok & McGaugh 1997 in LSB galax;- 
ies, Hernandez & Gilmore 1998 and Salucci 2001 for normal 
spirals, Tyson, Kochanski & DeU'Antonio 1998 for galaxy 
clusters). According to the analysis of rotation curves of 
dwarf and LSB galaxies carried out in Firmani et al. (2001), 



the core radius of the "virgin" halo of a galaxy of the size 
of the MW should be ~ 10 kpc with a central density of 
w 0.01 - 0.02 Mopc"^. To illustrate the effect that such a 
constant density core would have in terms of the rotation 
curve of the MW, we artificially smooth the inner density 
profile of the MW dark halo model along its evolution in 
such a way that at the present epoch it has such a core. The 
final rotation curve decomposition obtained after the formar 
tion of a MW disc in this halo is shown in Fig. 2a by the thin 
curves, which arc analogous to the thick curves of the NFW 
case. The rotation curve is now almost completely flat and 
the disc component dominates over the halo component up 
to the solar radius, i.e., the dashed and dotted curves now 
intersect at approximately the solar radius, as direct stud- 
ies of our galaxy show (e.g. Kuijken & Gilmore 1989). We 
can hence conclude that a comparison of dynamical stud- 
ies in the MW and galactic evolutionary models offers yet 
another evidence of the excessive central concentration of 
CDM haloes predicted by the standard cosmology. Fortu- 
nately, the inclusion of a shallow core in the halo has no 
major consequences upon the SFH, as will be shown below. 



4 STAR FORMATION HISTORY IN THE 
SOLAR NEIGHBOURHOOD 

As it is seen in Table 1, the properties of the MW fiducial 
model predict reasonably well the data inferred from obser- 
vations, with the exception of the rotation curve. In Fig. 3a 
we show the evolution of the SFR per unit area and of the 
gas infall rate, both at the radius iJo = 8.5 kpc, Ea(-Ro,i) 
(solid line) and Eg(_Ro,i) (point-dashed line), respectively. 
The disc at the radius Ro begins to form at a look-back time 
~ 11 Gyr {z « 2). 

An interesting test derives from taking into account an 
integrated colour index for the solar neighbourhood. Using 
recent Padova stellar evolutionary models (e.g. Girardi et 
al. 1996) we computed synthetic colour magnitude diagrams 
for the average SFH, shown by the solid curve in Fig. 3a, 
taking a solar metallicity for the last 2 Gyr, and one third 
solar before this point, as a first approximation to the en- 
richment history. Calculating the B — K colour index for 
realizations containing upwards of 200,000 stars we obtain a 
value of B — K — 3.15, in excellent agreement with the so- 
lar neighbourhood observational estimate for this quantity 
of B-K = 3.13 (Binney & Merrifield 1999). Repeating the 
experiment for SFHs resulting from MAHs deviating from 
the average one, we obtain practically identical values for 
the 7^— band luminosity, reflecting equal integrals under the 
SFH curves, but rather different values for the B— band lu- 
minosity, reflecting different recent SFHs. Taking the MAHs 
which differ most from the average (of the ones given in Fig. 
lb), we obtain differences of 0.3 mag in the predicted B — K 
of the solar neighbourhood, rather larger than the observa- 
tional uncertainties of the measured value. This important 
result leads us to expect the average MAH as the optimal 
choice for the Galaxy. 

The recent observational inferences of Binney, Dehen 
& Bertelli (2000 henceforth DBDOO) use kinematical data 
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Figure 4. Evolution of the relation between surface SFR and 
surface gas densities, at the given redshifts, for the solar neigh- 
bourhood and the integrated Galaxy, jointed triangles and circles, 
respectively. A Schmidt law of power 2 is clearly seen as the 
final result of the SF physics we used. The solid dots show the 
data of Kennicutt (1998), corrected to take helium into account, 
for which a Schmidt law of power ~ 2 is a good description. 

from the Hipparcos catalogue to break the age metallic- 
ity degeneracy at the oldest turn off, identified also from 
the Hipparcos satellite by combining photometric data with 
theoretical isochrones. Their results for the age of the so- 
lar neighbourhood are shown by the shaded bar in Fig. 3a, 
and are clearly in good agreement with our predictions for 
the average MAH. The maximum rate of gas accretion at 
Rq is attained at ~ 10.5 Gyrs, and strongly decreases after 
this epoch. The SFH attains a broad maximum at 8-6 Gyrs 
{z = 0.9 — 0.6) and then smoothly declines by a factor of 
« 2 towards the present day. The level of 3.1 Mopc"'^ Gyr^^ 
we obtain for this quantity lies well within the range deter- 
mined obscrvationally, shown by the vertical shaded box at 
t=0 Gyr (Talbot 1980, Rana 1987). A similar SFH is ob- 
tained for the MW model with a shallow core (thick dashed 
line), which illustrates the robustness of our results regard- 
ing SFHs, to the existing discrepancies in dark halo density 
profiles between theory and observations. 

Our SFR at the solax radius is almost indistinguishable 
from an average of an annulus with inner and outer radii of 
7.0 and 10 kpc, respectively. The local SFH inferred from 
observations might refer to an anrmlus of this size because 
due to stellar diffusion, stars which arc today in the solar 
neighbourhood could have been formed at Galactocentric 
radii differing by up to ~ 1.5 kpc (e.g. Dehnen & Binney 
1998b). 

The difference in shape between the SFH, Es(i?o, t), and 
the gas infall history, T,g{Ro,t), depends on the turbulent 
dissipation time ta; the shorter ta is, the closer T,siRo,t) is 
to Eg (Tio , t ) • Nevertheless, the dissipation time we use is self- 
consistently calculated according to Firmani et al. (1996), 



and is in rough agreement with results from MHD rmmerical 
simulations of the ISM (Avila-Reese & Vazquez- Semadeni 
2001). We note that matching the constraints on the age of 
the solar neighbourhood of BDBOO with a different value of 
td would require the use of a MAH differing largely from the 
average one, a rather improbable situation. 

In Fig. 3a we also show Eg for the solar neighbourhood, 
as inferred recently by Bertelli & Nasi (2001), jointed cir- 
cles. These authors used a volume-limited sample of field 
stars from the Hipparcos catalog and compared to synthetic 
color-magnitude diagrams generated for different parame- 
ter values of a SFH having a fixed parametric shape, and 
different stellar initial mass functions. Applying a mini- 
mization, the best parameters of the a priori given SFH are 
found. In order to pass from the SFR per unit volume to the 
SFR per unit area, the total number of stars of age t found 
within the sampled sphere, as a fraction of the total over 
the disc height, a time dependent correction factor, must be 
known. This can be obtained with a procedure which takes 
into account the details of the vertical disc force law and the 
variations of velocity dispersion with age (e.g., Hernandez et 
al. 2000b). 

The SFH shown in Fig. 3a corresponds to the solution 
for the model with more degrees of freedom of the ones used 
in Bertelli & Nasi (2001) (their var-var model with h = 2), 
and including the volume to surface density correction men- 
tioned above. It must be noted that the detailed shape of 
the observed SFH is the result of the age dependent correc- 
tion factor, as well as of the volume SFH solved for. As the 
authors emphasize, their method is useful to describe only 
the general trend of the SFH over the total life of the system, 
and not its detailed shape. 

As it can be concluded from the plots shown by these 
authors it is that the average SFR over the period 10-6 Gyr 
(in look-back time) could not have been larger than that 
over the 6-0 Gyr period. From Fig. 3a, this general trend is 
seen to agree well with the prediction of our fiducial model 
for the MW, suggesting that the main ingredients of the SF 
process in the MW were correctly talcen into account, and 
that the MAH of the MW has been close to the average one. 

The integration of the SFR in time gives the stellar 
surface density Eg, shown by the solid line in Fig. 3b, as a 
function of time, at Ro- The dashed line is the gas surface 
density Eg. At the present epoch, Ss(-R()) = 41.3 Mopc"'^ 
and Eg(_Ro) = 13.6 Mqpc^"^ in good agreement with ob- 
servational measurements (Table 1). The evolution of the 
stellar surface density at the solar neighbourhood obtained 
with our model compares well to that derived by Fuchs et 
al. (2000) from a sample of G and K dwarfs with individu- 
ally determined ages using their chromospheric HK emission 
fiuxes, and assuming a given gas accretion history (see their 
Fig. 7). 

Fig. 4 presents a global check of our SF mechanism. The 
open triangles show the time evolution of the solar neigh- 
bourhood in a log SFR vs. log Eg plot, at z=0, 0.5, 1 and 1.5. 
The open circles show the temporal evolution for the whole 
modeled galaxy at z=0, 0.5, 1, 1.5, 2 and 2.5. The solid 
dots show the observations of normal galaxies by Kennicutt 
(1998) corrected to take helium into account. We have not 
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Figure 5. a) SHFs resulting from the five random MAHs of Fig. 
lb, in matching hne styles. The thick dashed curve shows the 
indicative results of Bertelli & Nasi (2001). b) Here we illustrate 
the variations in the predicted SFH of the solar neighbourhood 
as resulting from changing the input value of A in consistency 
with the uncertainties in the observational determinations of the 
disc scale length, no strong dependence is seen. The thick dashed 
curve is the same as in panel a). 



included here the ULIRGs of that sample because they are 
not representative of the SF regime of normal galaxies but 
are related to very strong SF bursts. In these cases, when 
the SFR exceeds the critical value Es,c = 3OOM0pc~^Gyr~^ 
(Eg,c = IOOMqpc"^) then SN energy injection no longer 
controls the SF, the control is transferred to the radia- 
tion pressure (Eddington limit). This duality in the SF self- 
regulation mechanism has been studied by Firmani & Tu- 
tukov (1994). For normal galaxies where the SFR is typically 
less than Ss,c (the case of this work) self-regulation by SN 
energy injection leads to a Schmidt law of power ~ 2, while 
for the bursting SF regime with SFR greater than Ea,c the 
self-regulation by radiation pressure leads to a Schmidt law 
of power ~ 1. The comparison of the open and filled sym- 
bols shows a rather satisfactory agreement indicative of a 
Schmidt law of power 2. The slight offset of the data 
towards higher values of the SFR is not particularly mean- 
ingful because of the scatter of the observational data and 
the uncertainty on the turbulent dissipation time. 

The shape of our resulting SFHs for the whole MW disc 



differ from that at the solar radius. The inside-out formation 
of the disc and the mechanism of SF used here indeed lead 
to a differential evolution with radius of the galaxy (see also 
Boissier & Prantzos 1999, hereafter BP99), which highlights 
the importance of a dynamical galactic model in studies of 
galactic evolution where comparisons are made with spa- 
tially localized data. 

4.1 Sensitivity of the local SFH to the model 
input parameters 

The input parameters of the MW evolutionary model pre- 
sented above were chosen in such a way that the present-day 
properties obtained by the model agree with the observed 
structural parameters of the Galaxy, V50 , rg and Ms (see Ta- 
ble 1). Since the observational data are uncertain to some 
level and the model presents some degeneracy with respect 
to the input parameters, we need to explore the sensitivity of 
the predicted SFH in the solar neighbourhood to this input. 

Due to the difficulty in determining the MAH for the 
MW from any observational constrain (see §3), we have 
taken the average MAH (fiducial model) as a first choice, 
as it represents the most probable cosmological solution for 
a MW mass galaxy. This obviously presents a smooth shape 
(Fig. 2a). However, other MAHs are also able to produce 
dark haloes consistent with the MW. In Fig. lb we have five 
MAHs taken at random which, for a proper choice of the in- 
put parameters, can also give models in agreement with the 
main conditions we fixed for the MW. These MAHs are in- 
trinsically discontinuous due to the stochastic nature of the 
primordial density fluctuation field. The SFHs at Ro = 8.5 
kpc for the MW models produced with these MAHs are plot- 
ted in Fig. 5a, with the line style corresponding to the MAHs 
of Fig. 2a. We have included a smoothing on the results of 
these SFHs to facilitate the comparison with the results of 
Bcrtclli & Nasi (2001) (thick curve), as those results have a 
low time resolution. 

As it can be seen, both the age of the solar neighbour- 
hood and its present-day SFR arc sensitive to the assumed 
MAH of the Galaxy. It is therefore these dependences that 
allows us to constrain the MAH of the MW. The observed 
age of the solar neighbourhood and its present SFR imply 
that the MAH of our Galaxy has been very close to the 
average one for present day objects of that mass. 

The spin parameter A is mainly constrained by the disc 
scale length rg, which has an observational uncertainty of 1 
kpc. In Fig. 5b we present the SFHs corresponding to MW 
models with the average MAH but with rs=2.5 kpc (dotted 
lino) and rs=3.5 kpc (dashed line), the two extremes of the 
observational range. The difference in the SFHs is small. 
The more concentrated disc (smaller A) results in a later 
formation of the solar neighbourhood than the less concen- 
trated one (larger A). The maximum in the SFH thus shifts 
slightly in time for these two alternative models, although 
much less than the changes resulting from varying the MAH. 
Therefore, uncertainties in the scale length do not alter our 
conclusions regarding the connection between the MAH and 
the SFH of the Galaxy. 

Changes in fa and Mv allowed by the uncertainties in 
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V50 and Ms basically produce upward/downward shifts of 
the local SFH directly proportional to Mg, the shape of the 
SFH remaining almost the same. No large variations of these 
parameters are possible, if we want the present day Es and 
Eg at the solar neighbourhood to remain within the observed 
ranges. 

Finally, we note that the Hubble paxameter Ho plays a 
role in determining the cosmic time when the disc starts to 
form at the solar neighbourhood. The predicted SFH in the 
solar neighbourhood starts at lower look-back time when Hq 
is larger. For example, if we set Ho = 80km s~^Mpc~^, then 
the disc at the solar neighbourhood forms at « 9.0 Gyr, for 
the average MAH. 

We find that the disc age at Ro for the fiducial model 
with Ho = 65km s~^Mpc~^ is consistent with observational 
estimates of the age of the solar neighbourhood (BDBOO). 
Thus, the data offer a joint consistency check on the cos- 
mological scenario, the galactic evolution model and the SF 
physics used. 



5 DISCUSSION 

5.1 Comparison with chemo-spectrophotometric 
models 

In our approach, the local and global SF and infall gas his- 
tories of galaxies are predicted in a deductive way. In the 
case of models of chemical and spectrophotometric evolu- 
tion, these histories arc given in a parametric form and then 
constrained by comparing the model results with observa- 
tions. 

Much as the inversion of the high quality Hipparcos 
data allows for a reconstruction of the SFH of the so- 
lar neighbourhood, the large amount of observational con- 
straints both for the solar neighborhood and over the whole 
disc allow the cliemo-spcctrophotomctric models to constrict 
these histories (e.g., Prantzos & Aubert 1995; Chiappini et 
al. 1997; Prantzos & Silk 1998; BP99). However, most of 
the restrictions being time-integral constraints, the tempo- 
ral resolution attained is low. 

The SFHs in the solar neighborhood reported in Prant- 
zos &: Silk (1998) and in BP99 declines smoothly to the 
present epoch by a factor around two since its maximum, 
which is attained ~ 8 — 9 Gyr and ~ 6 — 7 Gyr ago, respec- 
tively. The main constraints for this smooth behaviour of 
the SFH, which agrees qualitatively with our prediction, are 
the observed differential metallicity distribution of long-lived 
G-dwarfs and the age-metallicity relation. These constraints 
also impose the necessity of a gradual gas infall, which again 
is in qualitative agreement with predictions of our model. 

However, due to the phenomenological nature of the 
chemo-spectrophotometric models, a detailed reconstruction 
of the SFH is difficult to attain, in particular at early epochs. 
For example, in both papers, Prantzos & Silk (1998) and 
BP99, the solar neighborhood is assumed to form 13.5 Gyr 
ago, which is too early according to accurate kinematical 
and photometric studies (BDBOO). Nevertheless, with am ap- 
propriate combination of the a priori gas infall history and 
SF law, the corresponding metallicity and surface density 



constraints have been roughly satisfied in spite of the early 
formation of the solar neighborhood (degeneracy). 

In as much as our results for global and local SFHs and 
gas infall rates do not deviate from those inferred by requir- 
ing consistency with the observational restrictions described 
above (BP99), we can expect our models to be also in agree- 
ment with the large variety of data regarding metallicities 
and abundance gradients in our Galaxy and for the solar 
neighbourhood. 

5.2 Stationary and bursting star formation 

Although wc have assumed throughout that the main trends 
of the SF processes in the disc are well described by the sta- 
tionary approach adopted, we are aware of the existence of 
ample evidence in favour of a more bursting character for 
this crucial ingredient, under some circumstances. However, 
as will be detailed below, we have reasons to believe that 
this alternative channel of SF will only introduce fluctua- 
tions of short temporal and spatial character, such that our 
stationary solution still represents the long term/large spa- 
tial trends accurately. 

The method applied by Bertelli & Nasi (2001) allows a 
reconstruction of the broad trends in the local SFH over a 
long time interval (~ 10 Gyr), but it is not useful to trace 
fluctuations of short duration. The question of significant 
temporal fluctuations has been raised by Hernandez et al. 
(2000b) who in a direct and statistically grounded study of 
the colour magnitude diagram of Hipparcos find a modu- 
lation of a factor of ^ 1.5, with a period of ~ 0.5 Gyr in 
the SFH of the solar neighbourhood. These authors inter- 
pret this as evidence of an enhancement of the SF activity 
caused by encounters with the spiral arm pattern density 
waves, a physical feature which is absent from our models. 
It is precisely the slowly oscillating character of the SFH 
found in that work, that justifies the interpretation of our 
results as the average over the inherent fluctuations of the 
Galactic SF process, i.e. this defines the temporal resolution 
of our present study to be of the order of 0.5 Gyr 

Even stronger time fluctuations in the SFH of the so- 
lar neighbourhood have been suggested by Rocha-Pinto et 
al. (2000), who argue against any steady (or self-regulated) 
component of the local SF. These authors use an empirical 
chromospheric activity indices vs. age relations to affix an 
age to each star in a sample of 552 late type dwarfs. How- 
ever, the empirical relation used by these authors presents 
a strong scatter of around a factor 2 in age (their Fig. 22), 
which should be included as a time smoothing kernel on their 
final results. Interestingly, if such a smoothing is applied on 
their SFH, a trend similar to that reported in Bertelli & 
Nasi (2001) appears (see Fig. 8 in Rocha-Pinto et al. 2000). 
We therefore find the above study in fact strengthens any 
conclusions drawn from the work of Bertelli & Nasi (2001), 
from which no objections to the SF physics we are using are 
apparent. 

A further mechanism responsible for strong time vari- 
ations in the SF processes of galaxies is that of merger 
and tidally induced SF enhancements. However, this applies 
mostly to systems such as the strongly interacting galax:ies 
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imaged by the IRAS satellite (Firmani & Tutukov 1994), 
or galaxies in dense cluster environments, see for a review 
Dultzin-Hacyan (1997). It is clear from the late type and 
low density environment of the MW that any alterations of 
the global SF due to interactions have been minor, and in 
all probability do not affect the conclusions of our study. 



6 SUMMARY AND CONCLUSIONS 

We have presented predictions of the SFH in the solar neigh- 
bourhood making use of a galactic evolutionary approach, 
tightly linked to a cosmological scenario. In this way we es- 
tablish a direct connection between the evolution of the so- 
lar neighbourhood and the cosmological background. Within 
our deductive approach, galactic discs form inside out within 
growing CDM haloes. Star formation is triggered by disc 
gravitational instabilities and is self regulated by an en- 
ergy balance in the turbulent ISM. In previous papers we 
have shown the ability of this approach to predict expo- 
nential discs over several disc scale lengths, as well as to 
reproduce the TuUy-Fisher relation and several of the cor- 
relations across the disc Hubble sequence. Our results are 
also in qualitative agreement with the gas and SFR histo- 
ries found to be consistent by comparing observations with 
chemo-spectrophotometric models. 

The allowed MAHs for the MW resuh in SFHs for the 
solar neighbourhood which bracket the range of observa- 
tional inferences. The average MAH for objects having the 
present day mass of the MW represents a cosmological max- 
imum likelihood solution for the Galaxy. The SFH in the so- 
lar neighbourhood predicted for the above possibility, in the 
flat A CDM model used here, closely matches the available 
observational data. 

For the cosmological galaxy evolution approach pre- 
sented, the age of the solar neighbourhood as inferred from 
observations implies a delay in the conversion of the accreted 
gas into stars, consistent with our physical formulation of the 
SF process in disc galaxies. 

For the value of Ho — 65km s~^ Mpc~^ used, the most 
probable MAH for our Galaxy implies an age for the so- 
lar neighbourhood of ~ 11 Gyr, in agreement with direct 
observational inferences. 

The rotation curve decomposition of the MW model is 
unrealistic if the disc forms within a CDM halo. If the halo 
has a shallow core of size and central density as inferred from 
rotation curves of dwarf and LSB galaxies, this problem is 
solved. Since the SFH is not affected by this, our conclusions 
regarding it and the cosmological build up processes of the 
Galaxy will remain for models where dark haloes present 
shallow cores. 
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